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1.  INTRODUCTION 

In  this  paper  we  shall  continue  our  study  of  various  sample  path 
properties  of  x2  processes,  begun  in  Aronowich  and  Adler  (1985). 

The  x2  process  is  defined,  in  law,  via  a  vector  X(t)  =  (X^(t), ...,X(t)) 
of  n  independent,  stationary,  zero  mean  Gaussian  processes  as 

n  ? 

(1.1)  Y(t) :  =  X  (t)X(t)  =  z  x:(t)  . 

i-1  1 

We  have  already  noted  in  the  previous  paper  why  this  process  is 
interesting,  and  we  refer  the  reader  to  the  introduction  there  for 
details.  For  the  sake  of  completeness,  however,  let  us  recall  here 
that  x2  processes  are  of  applied  interest  in  a  number  of  stochastic 
modelling  situations,  and  of  considerable  theoretical  interest  as 
perhaps  the  only  non-Gaussian  process  for  which  it  is  possible  to  develop 
a  theory  almost  as  general  as  in  the  Gaussian  case.  Thus  the  x2  process 
presents  a  natural  candidate  for  studying  the  robustness  of  the  Gaussian 
theory  when  similar,  but  non-Gauss ian ,  processes  are  considered.  As 
was  the  case  in  the  previous  paper,  our  principle  concern  will  be  to 
see  when  non-Gaussian  behaviour  is  exhibited  by  x2  processes,  rather 
than  trying  to  see  what  has  to  be  done  to  obtain  almost-Gauss ian  pheno¬ 
mena.  In  this  view  we  shall  pay  particular  emphasis  to  low  sections 
of  the  x2  processes  (which  is  clearly  bounded  below  by  zero)  which 
exhibit  highly  non-Gaussian  behaviour,  rather  than  emphasise  the  high 
sections,  which  are  much  closer  to  Gaussian.  This  emphasis  is  also  re¬ 
lated  to  and  motivated  by  the  applications  mentioned  in  the  earlier  paper. 
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The  paper  is  organised  as  follows.  Section  2  studies  the  fine 

behaviour  of  x2  processes  after  level  cross i ms  and  extrema.  Our 

approach  is  to  develop  the  Slepian  model  process  in  each  case,  this 
representing  the  conditional  distribution  of  the  original  process  after 
a  level  crossing  or  extremum.  The  model  processes  take  on  particularly 
simple  forms  if  the  level  crossing  or  extremum  is  either  at  a  high 
level  or  close  to  zero,  and  so  we  study  these  cases  in  detail.  In  all 

cases  it  turns  out  that  at  high  and  low  extrema  x2  processes  look 

like  random  parabolas,  as  do  their  Gaussian  counterparts  and  as,  in 
fact,  does  any  smooth  enough  function.  The  differences  between  the  y2 
and  Gaussian  cases  show  up  in  the  parameters  of  these  parabolas, 
which  are  generally  random  variables. 

From  the  simple  structure  of  these  limiting  parabolas  it  is  a 
simple  matter  to  deduce  the  distribution  of  variables  such  as  the 
duration  of  high  level  excursions.  Examples  of  this  type  of  calculation 
are  given  in  Section  3. 

In  Section  4  we  consider  the  asymptotic  Poisson  form  of  the 
(normalized)  process  of  downcrossings  of  a  low  level  by  a  x2  process. 
Upcrossings  have  already  been  considered  by  Sharpe  (1978)  and  Lindgren 
(1980a, b),  and  these  follow  the  Gaussian  model  quite  closely.  Down- 
crossings  of  low  levels  exhibit  somewhat  different  behaviour,  primarily 
because  of  the  fact  that  the  process  is  bounded  from  below.  As  is  well 
known,  results  of  the  types  of  Sections  2  and  4  are  often  combined 
in  practice  to  model  processes  as  a  sequence  of  random  parabolas  posi¬ 
tioned  over  the  points  of  a  Poisson  process.  Consequently,  the  results 
of  Section  4  are  of  reasonable  applied  interest. 
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Given  the  asymptotic  Poisson  result  it  is  a  comparatively  simple 
exercise  to  obtain  the  (normalized)  distribution  of  inf(Y(t):  te[u,T]> 
as  T  ~  .  This  is  done  in  the  final  section  5. 

We  conclude  this  section  with  some  definitions  and  notation. 

Throughout  the  paper  we  shall  assume  that  the  component  Gaussian 
processes  in  (1.1)  are  independent  and  have  a  common  distribution,  with 
covariance  function 


R ( t )  =  Rx(t):  =  E{X(s)X(s+t) }. 


We  also  assume,  without  any  loss  of  generality,  that  R(0)  =  1;  i.e. 
the  X^  all  have  unit  variance.  We  shall  assume  that  each  X..  is  twice 
continuously  differentiable  with  probability  one.  (c.f.  Cramer  and 
Leadbetter  (1967)  and  Dudley  (1973)  for  sufficient  conditions  on  R  for 
this  to  be  true.)  Thus  we  can  define  the  second  and  fourth  spectral 
moments  of  the  X^  ,  respectively,  as 


(1.2) 


t=0 


E  {  X?  ( t ) ) 


v  = 


v  : 
x 


where  X  and  X  are  the  first  two  sample  path  derivatives  of  X. 
Furthermore,  we  then  have  the  following  behaviour  for  P,  and  its 
derivatives  near  the  origin: 


(1.3) 


R(t)  =  1 


2  4 

At  x  vt 
2  4! 


+ 


o(t4) 


as  t  -*•  0  , 
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(1.4)  | R(t)  +  xj  <  0(t2)  as  t  -  0  . 

It  is  immediate  from  the  definition  (1.1)  of  the  x2  process 
Y  that  it  inherits  from  the  X^  both  stationarity  and  twice  differen¬ 
tiable  sample  paths,  so  that  we  can  express  its  first  and  second  order 
derivatives  via  those  of  the  .  We  have 

(1.3)  Y(t):  =  —  Y(t)  =  2  Z  X- (t)X.(t)  =  2X'(t)  X(t)  , 

dt  i=l  1  1 

where  X(t)  =  (X.(t) , . . . ,Xn(t) ) .  Similarly,  writing  X(t)  for  the 
vector  of  second  order  derivatives  of  the  X^  »  we  have 

(1.4)  Y(t):-£.Y(t)  -  2  "  [xf (t)  ♦  X  (t)X  (t)] 

dr  i  =  l  1  1  1 

=  2 {X 1 (t)  X(t)  +  x'(t)  X(t) } . 

Furthermore,  Y  has  mean  n  and  covariance  function 

(1.5)  Ry(t):  =  E{Y(s)Y(s+t) }  =  2nR2(t) 

To  avoid  uninteresting  technicalities,  we  shall  assume  throughout 
the  paper  that  all  the  joint  (finite-dimensional)  distributions  of 
Y,?  and  Y  are  non-singular. 

As  we  have  already  noted,  considerable  emphasis  will  be  placed  on  study¬ 
ing  the  sample  path  properties  of  x2  processes  via  model  processes 
created  by  conditioning  the  original  process,  in  a  horizontal  window 
sense,  on  some  event  having  occurred.  We  shall  assume  throughout 
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that  the  reader  is  familiar  with  this  idea,  and  with  the  basic  motivation 
and  methodology  behind  it.  A  full  survey  of  the  model  process  approach 
is  given  in  Lindgren  (1984),  and  a  background  sufficient  for  our  pur¬ 
poses  can  be  found  in  Section  10.3  of  Leadbetter,  Lindgren  and  Rootzen 
(1983). 

In  order  for  the  model  process  approach  to  be  meaningful  it  is 
necessary  to  assume  that  Y  is  ergodic.  This  clearly  follows  from 
the  ergodicity  of  the  processes,  which,  in  turn,  follows  from 

(1 .6)  R(t)  -*-0  as  ( 1 1  -*■  ®  , 

a  condition  which  we  shall  henceforth  assume. 

Finally,  in  order  to  compare  results  in  the  x2  case  to  corres¬ 
ponding  results  for  Gaussian  processes,  we  shall  require  a  Gaussian 
process  "matched",  in  terms  of  its  second  order  properties,  to  our 
X2  process.  Such  a  process,  Z(t)  say,  will  be  defined  to  have  mean 
n  ,  variance  2n  ,  and  covariance  function  (1.5),  identical  to  that  of 
the  x2  process . 
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> 


BEHAVIOUR  OF 


PROCESSES  AFTER  A  LEVEL  CROSSING 


In  this  section  we  study  the  x2  process  Y(t),  t  >  0,  given 
that  at  t=0  the  process  crossed,  either  from  above,  or  from  below, 
a  given  level.  It  is  now  well  known  that  the  appropriate  framework 
for  such  a  study  is  via  the  notion  of  "horizontal  window"  (h.w.) 
conditioning  introduced  by  Kac  and  Slepian  (1959),  and  the  so-called 
Slepian  model  process.  For  details  of  the  latter  we  refer  the  reader 
to  Lindgren's  (1984)  recent  review. 

We  commence  by  considering  Y(t)  ,  conditioned,  in  the  above 
sense,  on  having  an  upcrossing  of  the  level  u  >  0  at  time  t  =  0. 

We  write  Yu(t)  for  the  conditional,  or  "model",  process.  Then, 
writing  k  for  equivalence  in  law,  we  have 


Theorem  2.1 

(i )  Following  an  upcrossing  r 

(2.1)  Y  (t)  =  I  {/uR(t)V.  -  R(t)W . //a  +  A-(t)}2  , 
u  i=l  1  11 

where  the  random  vectors  V  =  (V^,...,Vn)'  ,  W  =  (W^,...,Wn)‘  have 
the  joint  density 


(2.2)  p(v,w)  =  7r,s-n2'(ri+1)/2r(n/2)(y  ’w)+exp{-w'w/2}, 


with  respect  to  (surface  measure)  *  (Lebesque  measure)  on  the  cylinder 
(y’v  =  l)  x  Fn.  Here  (•)  =max{*,0},  and  '^(t)  are  independent,  zero  mean, 
Gaussian  nrocesses,  independent  of  V,W,  and  with  the  common  covariance  function 


(2.3)  r(s.t)  =  R(s-t)  -  R(s) R(t)  -  ft(s)ft(t)/x  . 


(ii)  The  characterise  c  function  of  the  joint  distribution 
of  Yu(t1)>.. .»Yu(tk)  is  given  by 


(2.4)  *(vr  • • • *v^) 


u/2R-n/2  »  “  ufl 

/  ||z|exp{-  po- -  i  czldcdz 


2/2  7i  Xu  0 


where 


A  =  det 


B  =  det 


1  -2i  v-j 

-2i V2  Ri  2 

"2lvkRlk 

2i^R 

'2lvlR21 

1  -2i  \>2  • .  . 

"2ivkR2k 

2ir  R 

"2ivlRkl 

-2iv2Rk2  ... 

!-2ivk 

2i  -IR 

-2i v  j  R-j 

-2i ‘ ‘ 

■2ivkRk 

1 

2i 

2  i  ^2  ^  •  •  • 

2ivA 

-2u 

l-2iv1 

-2  i  v  ^  2  •  •  • 

'2lvkRlk 

R1 

_2ivlR21 

1  ~2 1  v  ^  * 

-2n  vkR2k 

R2 

_2ivlRkl 

_2lv2Rk2  ••• 

l-2i,k 

Rk 

_"2ivlRl 

“2 1  •  •  • 

*2ivkRk 

1 

-2U 

1 


and 


*i  -  R<v-  v  =  =  ^ft  ] 


Proof.  Part  (i)  of  the  Theorem  is  given  in  Lindgren  (1983c,  1984). 

The  distribution  of  V  and  W  is  the  h.w.  conditional  distribution  o 
X(0)//u  and  X(0)//a  ,  given  an  upcrossing  of  a  level  u  at  t=0 . 
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Note  that  although  X(0)  and  X(0)  are  independent,  V  and  W  are 

not . 

(ii)  Under  the  conditions  given  in  the  Introduction,  we  have: 

(2.5)  p(y  1 Y (0)  =  0,  Y (0)  10),  =  y'1/  |z|p(y,u,z)dz, 

n.w.  Q 

where  p(y,u,z)  is  the  joint  density  of  Y  =  (Y(t^ ) , . . .  ,Y(tk) )  1 , 

Y(0)  and  Y(0),  and  u  is  the  expected  number  of  upcrossings  of  the 
level  u  by  Y(t)  in  unit  time.  Thus 

(2.6)  y  =  (u/2)(n“1)/2  e‘U/2  * 


(see  Sharpe  (1978)) . 

The  joint  characteristic  function  of  Y,  Y(0)  and  Y(0)  is  given 

by: 

*(y,v0,5)  =  E(exp{i[y'Y  +  vQY(0)  +  cY(0)]>) 

k 

=  E  e xp ( i [  I  v  X'(t  .)X(t .)  + v  X* (0)X(0)  +  2 SX* (0) X(0 ) ] } 
j=l  J  J  -  J  u~ 

k  ?  ? 

=  IE(exp  i  [  Yj-X^Uj)  +  vQX^(0)  +  2-X](0)X1(0)])}n 
=  {E(exp  i[(X‘,  X1(0),X1(0))N(X^,X1(0),X1(0))'])}n  , 
where  X^  =  (X^  (t-j) , . . . ,X]  (tk) ) '  and 


N:  = 


Vg 

K  0 


—  •-  -  •-  ■  --  -  -  --  -  -  -  - — — • — ■ — - — — ^ ^ 
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The  covariance  matrix  of  the  normal  variates  ,X^  (0)  ,X^ (0)  is 
(c.f.  (1.2)) 


1 

R12 

Rlk 

-R1 

"21 

1 

^k 

"2 

-*2 

Rkl 

Rk2 

1 

A 

R1 

“2  • 

••  \ 

1 

0 

-»1 

■ 

A 

0 

X 

Therefore  (see,  for  example,  Lukacs  and  Laha  (1964)), 

*(y,v0,£)  =  ( de t ( I -2 i  MN))'n/2  , 
and  a  simple  calculation  shows  that 
det ( I -2i  M N)  =  A  -  2ivQB  , 

where  A  and  B  are  given  after  (2.4). 

Using  the  inversion  theorem  for  Fourier  transforms,  we  now  have 

k  o  °°  k 

(2.8)  p(y,u,z)  =  ( 2  7t  )  /  ...  |^(v,vQ,c)exp{-i  [-j‘y+v0u+cz]}dc  n  dv.  . 

-  <*>  -CO  j=0 

=  (20"k_2  f’“  f  e-i ~ /  e"i?z  /  (A-2iv0B)'n/2e'iUv«>dv0dc  n  dv. 

-CO  -CO  -CO  -oo  j =  1 


To  evaluate  the  innermost  integral  we  note  that 
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(2.9) 


co 

/(e+aivQ) 

-oo 


“J0  dvr 


jO  ,  a 

L27run/2'1e"u3/r(n/2)  ,  a 


+1 

-1 


for  Re(e)  >  0,  (Gradshteyn  and  Ryzhik  (1980),  p.  318.), 

Thus 

/(A-2iv0B)-"/2e-1u''Oduo  =  (2B)"n^2  /”(  ±-- 

—00  —oo 

=  27run/2"1(2B)'n/2exp(-uA/2B)/r(n/2) 


(since  Re(A/B)  >  0,  or  else  the  density  (2.8)  would  be  degenerate, 
which  would  contradict  the  assumptions  made  in  the  previous  section.) 

Substituting  (2.8)  and  the  above  into  (2.5)  and  changing  the 
order  of  integration,  we  now  obtain: 


p(y|Y(0)  >  O.Y(O)  >  o)h M  ■ 


•  (2-rk  / 


/e 


-v'y 


eu/ZB-n/Z 


2/2ttAu  6 


/  | z |  /exp[-  -  i ?z]d£dz  f  n  dv. 

o  -»  ^  3  j=i  J 


Since  the  expression  in  braces  does  not  depend  on  .y  ,  it  equals 
the  characteristic  function  ^(v)  of  (2.4),  and  the  proof  is  complete. 


★ 

Remark .  Similarly,  if  we  now  write  Yu  to  denote  the  process  Y  condi¬ 
tioned  (in  the  horizontal  window  sense)  on  a  downcrossing  of  Y(t)  of 

★ 

the  level  u  at  t=0,  then  it  is  easy  to  check  that  Yu  has  a  repre¬ 
sentation  identical  to  that  of  Yu  ,  with  the  single  difference  that  in 
the  density  (2.2)  the  expression  (v'w)+  should  be  replaced  by  (-v'w) 
Furthermore,  in  the  character!- Stic  function  (2.4)  the  range  of  the 
z-integration  should  be  changed  to  the  negative  half-line. 
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The  density  (2.2)  turns  out  to  be  relatively  awkward  to  work 
with,  primarily  because  of  the  function  (•)*  appearing  there. 

Even  the  evaluation  of  the  mean  and  covariance  functions  of  the 
model  process,  is  a  difficult  task,  since  it  involves  computation 
of  the  second  and  fourth  order  moments  of  V  and  W.  However,  the 
joint  characteristic  function  is  more  amenable  for  evaluation  of  the 
moments  of  Yu( t ) ,  and  leads  to 

★ 

Corol lary2 .1  The  model  process  Y  ( t) ( Y  (t))  has  the  mean  function 


(2.10) 


m(t)  'J 
m*(t) 


=  n  +  (u-n)R2(t)  +  R2(t)/A  *  R(t)R(t)  *£^u7T 


and  the  covariance  function 


p(s,  t)|  p  o  o 

(2.11)  >  =  2(n-2u)fT(t)IT(s)  +  2nlT(t-s)  +  4(u-n)R(t)R(s)R(t-s) 

P*(s,t)J 

+  (2/ x)(2u-iru-2)R(t)R(s)R(t)R(s)  +  (4/ x)ft(t)R(s)R(t-s) 


-  (2/A2)R2(t)R2(s)  -  /27Tu/A(2R(t)R(t)R2(s)+2R(s)R(s)R2(t) 
+  (1 / A) R(t) R( t) R2 (s )  +  ( 1 / A) R(s) R(s )R2 (t) -2R(s) R(t) R( t-s ) 


-  2R(t)R(s)R(t-s)}  . 

Proof.  The  results  follow  via  appropriate  differentiation  of  (2.4) 
with  k=l,2.  Tedious  calculations  are  omitted. 


In  order  to  compare  the  last  result  with  the  one  obtained  in 
the  Gaussian  case,  we  use  the  matched  Gaussian  process  Z(t)  men¬ 
tioned  in  the  Introduction  with  mean  and  covariance  function  equal 
to  those  of  Y(t).  For  this  process,  the  means  of  the  model  pro¬ 
cesses  at  upcrossings  and  downcrossings  are  given  by 


m2(t)l 

m*(t)j 


=  n  +  (u-n)R^(t)  +  /2nir/x  R(t)R(t)  , 


To  see  how  these  differ  from  (2.10),  Figure  1  gives  the  means 

of  the  model  processes  of  a  x2  process  and  a  matched  Gaussian 

process  conditioned  on  a  downcrossing  of  the  level  \  at  t=0. 

t2 

(Covariance  function  R(t)  =  e“^  .)  Note  that  the  x2  process 

has  a  flatter,  and  higher,  minimum  than  the  matched  Gaussian  process. 
Both  return  reasonably  quickly  to  the  mean  level  of  4. 


Figure  1  near  here 


We  now  turn  to  the  study  of  the  model  processes  and  Y* 
at  high  and  low  levels,  respectively.  In  these  cases  the  model 
processes  take  on  particularly  simple  forms  which  are  useful  both  in 
the  modelling  of  x2  processes  at  extrema  and  for  prediction  pur¬ 
poses.  We  need  to  start  with  a  technical  result. 


Lemma  2 . 1  Let  X  =  (X-j , . . .  ,Xm) 1 »  Y  =  (Yj,...,Y  )'  be  independent  random 
vectors;  X'vN(0,a2I),  Y  *  N(0,b2I).  Denote  X  =  X'X,  Y  =  Y'Y,  Z  =  X'Y  . 
the  joint  density  of  X,Y,Z  is  given  by 


(2.12)  p(x,y,z)  =  C(xy  -  z2) ^m_3^/2exp{ -x/2a2  -  y/2b2}  ,  x  >  0,  y  >  0, 

2  n 

xy-z  >  0; 

c_1  =  2m/7  (ab)mr(m/2)r((m-l)/2)  . 

Proof.  The  conditional  ch.f.  of  Y,Z,  given  X  =  x,  is 
^Y  z  |  X=x^  V1  ,v2  ^  =  'i  +  i  v  2  x  -  ^ ) 

m  n 

=  n  E(exp{iv1YJ  +  i^xkYk» 
k-1 

2  ni  2 

=  exp{-iv2  x 'x/4v-| }  n  E(exp{iv1(Yk  +  Vk/2V  })  . 


For  every  k=l,...,m  the  r.v.  Yk  +  v£Xk/2Vl  has  a  N(v2xk/2v-j  ,b2) 

distribution.  Therefore  the  expectations  in  the  last  expression  yield 

2 

the  ch.f.  of  a  non-central  x-|  distribution; 


iv2x2 

E(exp{iv1(Yk+v2xk/2v1)2})  =  0-2iv1b2)'1/2exp{— ^ ? - }  . 

■  4v-j  (1  -2ib  ^  ) 


Consequently 


^Y,Zj  X=x^vl *v2 


,)  =  (1  -2i v-,b  ) 


2, -m/2  _v"b  x'x 


exp{- 


2  2 
"2r 


2(1  -2i 


This  being  dependent  on  x  only  via  x'x  =  x,  we  can  replace  the  condi¬ 
tioning  event  accordingly.  Using  the  inversion  theorem  one  obtains 
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.2  2 

b  xv. 


p(y,z|x)  =  (2 „r2  J(l-2iv1b2)"m/2e'iviy  Jexp  -{ — 


y  -  i  Z\>0 }dv„dv-, 
2(1  -2  i  v  -j  b  ^ )  2  2  1 


.2  ,0.2 


_  exp{-z  /2b  x}  n  1  „_-,.2  N - ( m-1 ) /2  ,  . ,  2.  > 

- — 73/2 —  /(1_2lb  V  exp{-i  (y-z  /x) 

b/x(27r)'3/i  -  1 


v-|  }dv ^ 


Applying  (2.9)  to  the  above  yields 


p<y,zix) = ’ y  * j2/x  (x  *  o>  • 

0  p 

Finally,  (2.12)  follows  on  observing  that  X/a  is  a  variate. 


We  are  now  in  a  position  to  prove  the  following  result,  which  hinges 

on  the  representation  (2.1)  for  Yu  and  the  corresponding  representation 

for  Y*  . 
u 


Theorem  2.2  Define  the  normalised  and  limit  model  processes 
Yu(t):  =  Yu(t//u)  -  u  jYjt):  =  -At2  +  2/A  Zt 

Y*(t) :  =  (l/u)Y*(/u  t)  ;  Y*(t) :  =  An2  -  2/a  Zt  +  1 


where  Z  =  V’W  and  v  =  w’w.  Then  Z  and  v  are,  respectively  Rayleigh 
2 

and  xn+1  variables  with  joint  density 


,  ..  zh-zZ)ln-3)/2^n 


0,  >"  >  z 


Furthermore 


(  i)  PUim  sup  |Y  (t)  -  Y  (t)|  =  0} 
u^>|t|<t0  u 

(ii)  P{*im  sup  | Y*(t)  -  Y*( t) (  =  0} 
u»0|t|<t0  u  0 

Proof.  (i)  3y  (1 .3} 

R(t/ /u)  =  1  -  (1+o(1)) 

ft(t//5)  (1-K>0 )) 

/u 

as  u  oo,  where  the  o-terms  tend  uniformly  to  0  for  t  in  any  bounded 
interval.  From  Theorem  2.1  we  have 

Yu(t)  =  ~xt2(l+0(l))  +  2/xtV'W( l+o(1)) 
n  2 

+  I  {^(Wu)  +  2/uIl  -  (1+0(1)) ]V.Ai(t//iJ) 

+  2 t(l+o(l))WiA.(t//u)}  . 

By  (1.4)  and  (2.3)  the  A..(t)  have  continuously  differentiable 
sample  paths  with  A. (0)  =  A  (0)  =  0  ,  and  thus  a. (t//u)  =  o(t//u)  as 
u  ■+  °°,  uniformly  for  bounded  t.  It  remains  therefore  to  find  the  dis¬ 
tribution  of  Z  =  V'W.  But  since  this  is  equivalent  to  the  h.w.  conditional 
distribution  of  Y(0)/2/uX  ,  given  an  upcrossing  of  a  level  u  at  t=9, 
its  density  is  given  by 

(2.13)  p(z)  =  vj_14uAzp(u,  2/u ~  z)  ,  z  >  0, 


=  1 

=  1  . 


where  p(-,0  is  the  joint  density  of  Y (0 )  and  Y(0)  .  Given 
X(0)  =  x,  ?(0)  =  2x'X(0)  is  a  normal  zero  mean  r.v.  with  variance 
4ax'x  ,  and  since  this  depends  on  x  only  via  x'x  =  u,  we  may  conclude 
that  Y(0) ,  given  Y(0)  =  u,  has  (conditionally)  a  N(0,4Au)  distribution. 
Therefore 


P(u,*)  =  p(y|u)p(u) 


u(n-3)/2e-u/2exp{_^2/8Au} 
Ax  2{n+3)/2T(n/2)  . 


and  a  direct  substitution  into  (2.13)  leads  to 


p(z)  =  z  e 


z  >  0  , 


the  desired  result. 

(ii)  In  a  similar  way,  via  (1 .3)  and  (1 .4),  as  u  ->  0 
R( AT  t)  =  1  -  £  xut2(l+o(l)) 

R(Ai  t)  =  -aAT  t(l+o(l)) 

(AT  t)  =  o ( /IT  t) 

uniformly  for  bounded  t.  From  the  remark  following  Theorem  2.1  we  can 
then  deduce  that 

Y*(t)  =  1  +  xt2W'W  +  2/AtV'W  +  o(/u)  . 


The  joint  distribution  of  Z  =  V'W  and  y  =  W'W 
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is  equivalent  to  the  h.w.  conditional  joint  distribution  of  ?(0)/2/uA 
•  • 

and  X'(0)X(0)/A  ,  given  a  downcrossing  of  the  level  u  at  t=0  .  Their 
joint  density  is 

(2.14)  p(z,^)  =  u  ^4uX | z | p(u»2>^uX  z,  \p)  ,  z  <  0 


where  p(-,-,-)  is  the  joint  density  of  Y(0),  Y(0)  and  X'(0)X(0)/x  . 
By  Lemma  2.1,  this  density  is 


p(y.y.'p) 


n-3 

(m  -  T /4)  exp{-y/2  -  »/ 2} 

2*+)A  An/2_1r(n/2)r((n-l)/2) 


and  a  substitution  into  (2.14)  concludes  the  proof. 


The  above  theorem  can  be  summarized  as  follows.  Near  an  upcrossing 

2 

of  a  high  level  u  the  x  process  can  be  approximated  by  a  random 
parabola 

-Aut2  +  2/Au  Zt  +  u  . 

2 

Near  a  downcrossing  of  a  low  level  u  the  x  process  can  be  approximated  by 

9  _ _ 

ah  -  mT zt  +  u 

which,  because  of  the  random  coefficient  ¥,  is  of  quite  a  different  form  to 
the  Gaussian  case. 

(It  is  interesting  to  note  that  the  model  process  in  the  Gaussian  case 
takes  a  similar  asymototic  form,  viz. 

-i-Aut2  +  AZt  +  u, 
with  Z  once  again  Rayleigh.) 
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It  is  intuitively  obvious  that  low  and  high  level  crossings  should  be 
followed  by  extrema.  The  model  processes  for  extrema  should  therefore  have  a 
very  similar  form  to  those  for  crossings.  We  shall  derive  the  1  imiting  model  pro 
cesses  of  extrema  with  the  aid  of  the  following  result. 

Lemma  2.2  Let  )(  =  (X(t^ ) , . . .  ,X(tj<)) 1  where  X(t)  is  a  stationary, 
Gaussian,  zero  mean  process  with  covariance  function  R(t)  (R(0)=1) 
and  finite  second  and  fourth  spectral  moments  (x  and  v  respectively). 


Then  the  conditional  distribution  of  X,  given  X(0)  =  x,  X(0)=x,  X(0)='x 
is  multivariate  normal  with  expectation 


.t  *■*?**!  ,  *t , .  V  y.l 

i  ji=l . I 


and  covariance  matrix 


(p(t.,t.)[  =  )r..  -  rr-TT  [AaR.R.»R.R.+R.R. 

I  1  J  )i,j-l,...#k  I  1J  1  J  1  J  J  1 


+  (a-l)R. R-  +  —  r.  r .] 
v  I  J  X  1  Jj 


i » j  1 , . . . ,  k 


where  a  =  v/x  . 


Proof.  The  covariance  matrix  of  X,X(0) ,X(0) ,X(0)  is 
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From  the  basic  properties  of  the  multinormal  distribution  it 
follows  that  the  conditional  distribution  of  X  is  normal  with  mean 

Y  ~  S .j 2^22 ^  ^ ’  ^0 

and  covariance  matrix 

{p(ti»tj)}i J  =  S11  "  S12S22S21  * 

Remark.  Lindgren  (1970)  derived  a  similar  result  for  x  =  0  . 

Applying  the  arguments  of  Lindgren  (1984)  it  follows  from  the 
above  lemma  that  the  model  process  obtained  by  conditioning  Y(t) 
on  having  a  local  extremum  of  height  u  at  t=0  can  be  written  as 

n  ..  9 

(2.15)  E  {/u  R(t)V.  -  R(t)W.//x  +  [ XR(t)  +  R(t)]H./x4T  +  K.(t)T  , 
i=l  1  1  11 

where  (t)  are  independent,  normal,  zero  mean  processes,  independent 
of  V,  W,  H,  and  with  the  common  covariance  function  p(s,t),  of 
Lemma  2.2.  The  distribution  of  V,  W,  H  is  the  h.w.  conditional 
distribution  of  X(0)//u  ,  X(0)//x  and  [X( 0 )  +  AX(0)]/Av/a-l  , 
given  an  extremum  of  hei  ght  u  at  t=0  .  Note  that  with  probability  1 
V*  V  =  1  and  V'W  =  0  .  Also,  /u( a-1 )  V'H  +  W'W  -  u  is  restricted 
to  negative  or  positive  values,  respecti vely, according  to  whether  (2.15) 
comes  from  conditioning  on  a  maximum  or  minimum.  It  was  shown  in  Arono- 
wich  and  Adler  (1985)  that,  given  Y(0)  =  u  and  Y(0)  =  0  ,  the  random 
variables  X ' (0) X (0 ) / x  and  X'  (0)X(0)/  -Vu( V-l)  are  conditionally  inde¬ 
pendent  having  x?  i  and  ~1~~T  +  N(0,1)  distrinutions  respecti vely . 


Furthermore,  since 


R(t)  =  -A  +  (v/2)t2  +  o(t2)  ,  t  -*  0 

and 

K.(t)  =  o(t),  t  -  0 

with  probability  1,  we  have  the  following 

■flp 

Theorem  2.2  Let  y  (t)  and  ^y(t)  be  the  h.w.  model  processes  giver 
maxima  and  minima  of  height  u  ajt  t=0  respectively.  Write 

yu(t)  =  yu(t//0)  -  u 

y*(t)  =  (l/u)/*(^7  t)  . 

Then 

(  i)  PUim  sup  I  V  (t)  +  At2j  =  0}  =  1 
u-*»|t|<t0  u 

(ii)  P{tim  sup  I  y*(t)  -  Ayt2  -  1|  =  0}  =  1  , 
u+0 j  t|  <tQ 

where  y  =  W'W  is  a  x2n+^  random  variable. 

Proof.  Using  the  preceding  remarks,  the  proof  is  strai ghtforward. 

3.  THE  limiting  distribution  of  length  and  height  of  an  excursion 

As  a  consequence  of  Theorem  2.1  we  can  obtain  the  limiting  di stri 

bution  of  the  length  and  height  of  an  excursion  above  high  and  below 

low  levels.  Let  t  and  t*  be  the  lengths  of  h.w.  excursions  by 

u  u  3  J 

Y(t)  above  and  below  a  level  u,  respectively,  i.e. 


tu  =  tnf{t  >  0:  Yu(t)  =  u} 
t*  =  inf{t  >  0:  Y*( t)  =  ul. 


Theorem  3.1 


(  i)  P{Him  /u  t  =  2Z//T}  =  1. 
u-*»  u 


( i  i )  PUim  x*//u  =  }  =  1  . 


Proof.  From  Theorem  2.1  we  have  that  with  probability  1 


«-im  /u  t  =  Him  inf{t  >  0:  Y  (t//u)=u}  =  2Z//x 
u-*x>  u  o+«>  u 

Him  t  J/u  =  Him  inf{t  >  0:  Y*( /lT t  )  =  u}  = 


By  using  the  joint  density  of  Z  and  v  ,  it  is  easy  to  obtain  the 

97 

density  of  5  =  as 


=  ,Fl(n±3  ,  n+l.  -  ^  )  , 


x(n*l)/2n,5n+2 


E  >  0 


where  ^ is  a  degenerate  hypergeometric  function. 

★ 

Now  let  Qu  and  be  the  heights  of  h.w.  excursions  by  Y(t) 
above  and  below  a  level  u,  respectively.  Thus 


Qu  =  rrax{Yu(t)  -  u:  0  <_  t  <_  tu) 


Q.  =  minfu  -  Y  (t) :  0  <  t  <  t  1  . 
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Theorem  3.2 

(  i )  P{  £i m  Q  =  Z2 }  =  1 . 

U-*x= 

(ii )  P{  £im  Q*/ u  =  }  =  1 . 

Proof.  We  have,  with  probability  1, 

tim  Q  =  Jlim  max  {Y  (t/Zu)  -  u:  0  <  t  <  2Z//x  }=  , 

U-x»  U-x»  U 

Him  Q  /u  =  Jlim  min{l  -  Y  (/u  t ) / u :  0  £  t  £  }  =  Z^/t 

u-t-0  u  u+0  u  /if 


2 

Note  that  Z  is  an  exponential  r.v.  with  mean  2,  and,  as  a  simple 

2 

calculation  shows  ,  x.  =  Z  h  has  the  density 


PU)  =  ^  ^n"3^2  ,  0  £  z,  <_  1 


4.  ASYMPTOTIC  POISSON  CHARACTER  OF  THE  PROCESS  OF  LOW  LEVEL  DOWN CROSSINGS 

It  is  well  known  that  the  high  level  upcrossings  and  (by  symmetry) 
the  low  level  downcrossings  of  a  Gaussian  process  often  behave,  after 
appropriate  normalization,  like  the  points  of  a  Poission  process.  This 
occurs  under  two  conditions.  One,  a  local  condition,  requires  that  the 
sample  paths  of  the  process  be  sufficiently  smooth.  The  second  is  a  mixing 
condition,  generally  expressed  by  requiring  that  the  covariance  function 
tends  to  zero  fast  enough  at  infinity. 

2 

Sharpe  (1978)  extended  the  upcrossing  result  to  >■  processes  under 
the  mixing  condition  R(t)  =  0(t  a)  for  some  0  <  a  <  1  as  t  *  *>,  and 
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Lindgren  (1980)  sharpened  this  by  requiring  only  ,R(t)nn  t  -*■  0  as 
t-*-®.  In  this  section  we  shall  look  at  the  corresponding  result  for 
downcrossings  of  low  levels.  (Unlike  the  Gaussian  case,  this  does  not 
follow  trivially  from  the  high  level  result.)  As  did  Sharpe,  we  shall 
model  our  proof  very  closely  on  Cramer  and  Leadbetter's  (1967)  presentation 
of  the  Gaussian  case.  In  the  end  we  shall  require  the  stronger  mixing 
condition  R(t)  =  0(t*a)  rather  than  the  weaker  R(t)  £nt  -+  0  . 

We  note  now  that  this  is  not  due  to  the  "older"  style  of  proof,  but  is  due 
to  the  fact  that  the  left  hand  tail  of  a  x2  distribution  has  a  power 
decay  to  zero,  unlike  the  exponential  decay  of  the  riqht  hand  tail.  We  also 
note  that  although  Berman  (1983,1984)  has  a  general  theory  of  rare  events 
that  would  seem  to  incorporate  the  next  result,  it  turns  out  that  verifying 
that  our  conditions  fit  into  his  framework  is  as  much  work  as  proving  every¬ 
thing  virtually  from  scratch. 

To  avoid  needless  repetition  of  considerable  detail,  we  shall 
assume  throughout  the  following  that  the  reader  is  equipped  with  copies 
of  both  Sharpe  (1978)  and  Cramer  and  Leadbetter  (1967). 

Theorem  4.1  Let  the  common  covariance  function  R(t)  of  the  component 
Gaussian  processes  X^(t)  ,  i=l,...,n,  sati sfy 

(4.1)  R(t)  =  1  -  2l  +  5T  t4  +  0U4)  as  t  •*  0, 

(4.2)  R(t)  =  0(t"a)  for  some  0  <  a  <  1  as_  t-+®. 

Let  Du(t)  denote  the  number  of  downcrossings  by  Y(t)  of  the  level  u 

in  the  interval  (0,t),  and  let  y  =  EDu(l)  .  Then  for  every  fixed  k=0,l,..., 

k  -t 

(4.3)  £im  PCD  (  L  )  =  k]  =  T-£r-  . 

u  u 
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Proof.  The  (Gaussian)  proof  in  Cramer  and  Leadbetter  (1967)  is  set  out 
in  the  form  of  five  lemmas,  three  of  which  carry  over  with  only  minor 
modifications  in  notation:  Lemma  1  uses  merely  the  fact  that  the  limit 
of 

u  -  F  (  u  i  e'u/2 

{4*4)  y"/7(2)  7(i72) 


is  zero  as  u  -*■  0;  Lemma  3  follows  directly  from  Lemma  2;  while  the 
proof  of  Lemma  5  is  a  direct  adaptation  of  Sharpe's  (1978)  proof  of  his 
result  (4.20) . 

We  are  left  therefore  with  the  task  of  proving  Lemma  2  and  pointing 
out  the  differences  in  the  proof  of  Lemma  4. 

Lemma  2  can  be  expressed  as 


(4.5)  Qu  =  E{Ou(t])CDu(t1)-l]>/EDu(t1) — -  0 

u^O 

for  an  interval  of  length 


t 


1 


-xLlIl. 

utx/u8*1 


] 


where 

(4.6)  0  <  (k+4)e  <  a  <  1. 

(For  notation  see  Cramer  and  Leadbetter  (1967)). 

Proof  of  (4.5) .  Using  (4.5)  of  Sharpe  (1978)  and  (12.2.10)  of  Cramer 
and  Leadbetter  (1967)  we  write 
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(4.7) 


1  1 

4-/  ds.  /  ds?E[r(s1)Y(s2)|Y(s1)=Y(s2)  =  u]p(u,u)  , 
0  0 


where  p( - , * )  is  the  joint  density  of  Y(s.|)  and  Y(s2)  .  The  density 
p(u,u)  can  be  obtained  by  a  strai ghtforward  modification  of  the  expres¬ 
sion  given  in  Miller  (1980),  p.  162  for  the  2-variate  Rayleigh  density 
function 


u|r| 


(4.8)  p(u,u)  = 


n_  _  i  _u _ 

u2  e  1-R2  In_2  (1-R^ 
~T~ 


- rT7~ 

(1-R2)  |  R|  2  r(§-)2? 


+ 1 


I 


where  R  stands  for  R(s^-s2)  and  Iv  is  the  modified  Bessel 
function  of  the  first  kind  and  order  v. 

Substituting  (4.8)  and  the  result  (Sharpe,  p.  379) 


e[y-(Si)Y'(s2)I'  vSiH(s2)-u]  i2u[a(1-R2)-R2]/(1-R2) 


into  (4.7),  we  find 
t 


(4.9)  1-R  >U-R  )„-R.  __  ^  . 

0  0  <i-rz)zi R|~r~  ~r~  l_R 


with  K  an  unspeci fied  posi ti ve  constant. 

Making  the  change  of  variables  s  =  s2~s1  and  denoting  the  inte¬ 
grand  in  (4.9)  by  g(s)  we  obtain 
tl 

Qm  1  K  /  g(s)ds 
u  0 

by  the  virtue  of  g  being  an  even  function. 
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An  asymptotic  expansion  of  I  for  large  arguments 
(Abramowitz  and  Stegun  (1964),  p.  377)  gives  that  for  some  x  >  1 
there  exists  a  C  >  0,  such  that 


I  (x)  <  Cex//x  . 

Taking  x  =  ,  we  have  that  as  long  as 

(4.10)  1  >  |  R(s)|  >  -  1  , 


then 


(4.11)  I  jM)  <  C(^-1/2exp{iil4 
^-1-r  l-RT  1  -FT 


Let 


Sq  *  inf  {s  >  0:  |  R(s)|  <  -  ^  +  1  } 


From  condition  (4.1)  we  deduce  that  there  is  a  positive  such  that 

(4.12)  1-C^2  >  | R(s)[  for  0  <  s  <  sQ  . 


Write 


;  *  - 


)] 


1/2 


Clearly  sQ  <_  s*  and  s*'Ik7  as  u  0. 


1  -  C}sZ/Z  <  (1  +  C^2^)"1  , 


Since 


it  follows  from  (4.12)  that 


(4.13)  1  +  |R|  <  2  -  C,s2  <  — £—3 —  for  0  <  s  <  sr 

1  l+C1$2/2  C 

Using  (4.1),  it  can  be  shown  that 

(4.14)  =  —  s  +  o(s)  as  s  •>  0  . 

(1-RZ)3/Z  2/x 

Applying  (4.11),  (4.13)  and  (4.14)  we  have 

S°  S°  l  J  ux-l  u 

/  g(s)ds  <  Kf  s(*4  +  t--t)  exp  -  {— - }ds 

0  0  4  1  +  |R| 

<K5i[»(s*/c^7r-  ) 

^  K[<j>(  u/2)  -  i-]  0  as  u  -*•  0  , 


where  $(•)  is  the  standard  normal  probability  integral. 

Furthermore,  since  exp{-  u  *  }  is  bounded  away  from  zero,  and 

1-R 

->  >  ass*.  , 

o-R2r 


(this  follows  from  (4.2)),  we  have  for  the  remaining  part  of  our 


>  * 


where  Z(s)  =  —  ^  ^  - 
1-R  (s) 


Since 


z  Iv(z)  -*-  2 /r( v+1 )  as  z  -*■  0  , 


(Abrahmowitz  and  Stegun  (1964,  p.  375),  the  last  inequality  becomes 
t1  n-1 

/  g(s)ds  <  KC^-SqJu"^ 
s0 

B-l 

From  t-j  ^  v  as  u  ->  0  we  have  for  sufficiently  small  u  that 
t1  <  2y6-1  ,  which  leads  to 

t  n_1 

|  ^(sjds  <  Ku~^~  (2ye_1  -  sn) 


0  as  u  -*  0. 


This  completes  the  proof  of  Lemma  2. 

The  only  modification  in  the  proof  of  Cramer  and  Leadbetter's  (1967) 
Lemma  4  to  be  made  is  the  order  of  the  probability  of  the  event  gr 
(pp.  265-266).  This  event  in  our  case  becomes 


9r  =  (Y(vrq)  <.  u} 


of  probability 


P(9r}  =  0(un/2)  . 


To  avoid  confusion  we  will  denote  Cramer  and  Leadbetter's  n  by  n' 
Consequently  the  final  part  of  Lemma  4  becomes 
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P(Dk>  -  P{Ek>  =  0(n 


,k+1un/2)  =  0(u2 


8< 


l-B(k+l))+B(k+l)/2 

) 


->  0  as  u  +  0  , 


where  of  course,  in  the  definitions  of  the  above  events,  upcrossings 
are  replaced  by  downcrossings  and  U(vq)  >  u}  by  (Y(vq)  <  u}  . 

5.  LIMITING  DISTRIBUTION  OF  EXTREME  VALUES 

The  limiting  extreme  value  distribution  of  max  Y(t)  as  T  -+  »  was 

0<t<J 

derived  by  Sharpe  (1978)  who  showed  that 

P[l-max  Y(t)  -  uT  <_  z]  -»•  e"e  as  T  +  ® 

Oft<J  1 

where  n_-j 

ut‘  lost^rTOogT)”]  . 

Thus  the  asymptotic  distribution  of  max  Y(t)  is  of  Type  1  (c.f.  Lead- 

0<t<T 

better,  Lindgren  and  Rootzen  { 1983) ), precisely  as  in  the  Gaussian  case. 

One  expects  that  since  x2  processes  are  bounded  below  ,  with  power-type 

densities,  the  asymptotic  distribution  of  min  Y(t)  should  belong  to  a  differ- 

0<t<T 

ent  domain  of  attraction.  This  is  in  fact  the  case,  as  the  following  result  indicates 


Theorem  5.1  The  minimum  of  a  x2  process  has  an  asymptotic  distribution 
of  Type  3.  Speci fi cal.y,  under  the  conditions  of  Theorem  4.1 
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where 

,  fTh  ~r 
(5.2)  VT  =  2  C  r(n/2) 

Proof.  It  follows  from  the  asymptotic  Poisson  character  of  the  down- 
crossings  process  that 

P{min  Y(t)  >  u}  ■>  e"T  as  u  -*•  0  , 

0<t<T 


where  u  and  T  are  related  via 

rrn 

<5-3> 

for  large  T  and  small  u  . 
n-1 

Putting  t  =  z  in  (4.3)  we  obtain 

u  2  z  , 

T/  A  / 7T 


and  (5.1)  -  (5.2)  readily  follow. 
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Mean  functions  for  >,^  and  matcnec  Gaussian  process? 
in  the  neighbourhood  of  a  downcrossing  of  H  at  t=C 


